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\ We present the results of two complementary investigations into the nature 

of strong (rest equivalent width, W r > 1.0 A) Mg II absorption systems at high 
redshift. The first line of questioning examines the complete Sloan Digital Sky 
Survey Data Release 3 set of quasar spectra to determine the evolution of the 
incidence of strong Mg II absorption. This search resulted in 7421 confirmed 
.£h ! Mg II systems of W r > 1.0 A yielding a > 95% complete statistical sample 

of 4835 absorbers (systems detected in S/N > 7 spectral regions) spanning a 
redshift range 0.35 < z < 2.3. The redshift evolution of the comoving line- 
of-sight number density, ^Mg(A) is characterized by a roughly constant value 
at z > 0.8 indicating the product of the number density and gas cross-section 
of halos hosting strong Mg II is unevolving at these redshifts. In contrast, one 
observes a decline in ^Mg(A) at z < 0.8 which we interpret as a decrease in the gas 
cross-section to strong Mg II absorption and therefore a decline in the physical 
processes relevant to strong Mg II absorption. Perhaps uncoincidentally, this 
evolution roughly tracks the global evolution of the star formation rate density. 
Dividing the systems in W r sub-samples, the ^MgPO curves show similar shape 
with lower normalization at higher W r values and a more pronounced decrease 
in ^MgPO at z < 0.8 for larger W r systems. We also present the results of a 
search for strong Mgll absorption in a set of 91 high resolution quasar spectra 
collected on the ESI and HIRES spectrographs. These data allow us to investigate 
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the kinematics of such systems at 0.8 < z < 2.7. In this search a total of 22 
systems of W r > 1.0 Awere discovered. These systems are characterized by 
the presence of numerous components spread over an average velocity width of 
Av ~ 200 km s _1 . Also, absorption due to more highly ionized species (e.g., 
Al III, C IV, Si IV) tend to display kinematic profiles similar to the corresponding 
Mg II and Fe II absorption. We consider all of these results in light of two 
competing theories previously introduced to explain strong Mgll absorption: 
post-starburst, supernovae-driven galactic winds and accreting gas in the halos 
of massive galaxies. The latter model is especially disfavored by the absence of 
evolution in £ Mg (X) at z > 1. We argue that the strong Mg II phenomenon 
primarily arises from feedback processes in relatively low mass galactic halos 
related to star formation. 

Subject headings: quasars: absorption lines (Mgll) 

1. Introduction 

Starbursting (SB) regions are likely to host rapidly expanding, highly ionized gaseous 
bubbles ("superwinds"). These bubbles are metal-enriched and powered by the correlated 
supernovae explosions of the stars on the high end of the SB region's mass-function. If the 
starburst is related to a merger event, shock- heating during the interaction may also lead to 
a wind of ionized gas (Cox et al. 2004). When nested in low enough mass galaxies, these 
bubbles have a good chance of overpowering their host's potential well and escaping to pollute 
the intergalactic medium (IGM) (Furlanetto & Loeb 2003; Madau et al. 2001), though 
Cox et al. (2004) note that even at high mass merging galaxies may have sufficient energy to 
produce expanding, gaseous bubbles. At high redshift, such bubbles are expected to travel 
at velocities of 100 to 1000 km s _1 and eventually expand to radii of ~ 100 kpc. Presumably 
the winds entrain the metals of prior and recent supernovae activity; therefore, the bubbles 
may cause metal-line absorption in background sources in the 0.5 < W r < 5 A range. 
Indeed, at low redshift the dwarf starbursting galaxy NGC 1569 is observed to possess these 
characteristics (Martin et al. 2002), as is NGC 1705, another dwarf starburster (Heckman et 
al. 2001). This enrichment process is believed to have played a leading role in the production 
of the observed metallicity of the present-day IGM (Madau et al. 2001; Scannapieco et al. 
2002), and therefore further observational understanding of this pollution mechanism is key 
to the refinement of enrichment scenarios. Furthermore, an exploration of the starbursting 
regions in absorption will constrain dynamic models of starburst winds and serve as an 
observational bridge for understanding the starbursting process at high redshift. 
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A link has recently been proposed between these metal-enriched expanding bubbles 
and the very strong Mg II absorption systems (W r > 1.8 A) observed in quasar spectra 
(Bond et al. (2001a), hereafter BCCV). In this paper the cutoff W r has been relaxed to 
> 1.0 A, as there is no known specific cutoff to the absorption depth of these systems and, 
as demonstrated later, smaller systems display features consistent with this phenomena. It 
is, however, demonstrated that stronger systems are more likely to be associated with this 
absorption host scenario. Assuming a simple spherical shape for the expanding shell, the 
expected profile for such a system is two roughly symmetric absorption lobes, corresponding 
to the cooler, expanding regions of the bubble, one in front and one behind the host galaxy 
along the line of sight, centered around a kinematically small, hot, non- absorbing region (a 
spectral "inversion"), where the actual host is located. As the expanding gas is subject to 
Raleigh- Taylor instability, it is expected to clump, causing distinct kinematic components to 
appear in the absorption spectra. Cloud-to-cloud metallicity variation may also be expected, 
a result of unique supernova history and/or incomplete mixing. The expected velocity range 
of such systems is ~ 100 — 1000 km s -1 . BCCV identified all of these features in the profiles 
of a small set of strong Mgll absorption systems. Similarly, Steidel et al. (2002) find that 
gas "infall" from such winds could help explain the kinematics they observed in the spectra 
of a strong Mg II system. These authors also emphasized that the study of more highly 
ionized species (e.g. CIV) may further reveal the origin of the Mgll gas and its relation to 
SN superwinds (see also Bond et al. 2001b). 

Based on the assumption that strong Mg II absorbers trace superwinds, and therefore 
starbursting galaxies, BCCV make a crude prediction of the redshift path density (dN/dz) 
of Mg II absorption. These authors found it to be consistent with a value extrapolated 
from the Mg II survey presented in Steidel and Sargent (1992), hereafter SS92, though the 
comparison suffers from small number statistics. Several goals of the work presented here 
include an extension of the work of BCCV to higher redshift, to obtain a large statistical 
sample, and to investigate evolution with redshift. 

An interesting correlation could also be made between strong Mgll absorption systems 
and Lyman Break Galaxies (LBGs), which are thought to frequently harbor starbursting 
regions. Comparisons of the velocities of the interstellar absorption features with nebular 
emission lines suggests wind speeds of 200 — 1100 km s" 1 (e.g. Pettini et al. 2001). In the 
gas, one identifies metal-line transitions of both low and high-ion species with significant 
column densities. Therefore, one may probe the physical characteristics of these winds (e.g. 
speed, cross-section, ionization state, etc.) through absorption- line studies. In turn, these 
observations should constrain the incidence and evolution of superwinds during the first 
~ 5 Gyr of the universe. 
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A competing scenario for the origin of strong Mgll absorption comes from Mo & 
Miralda-Escude (1996), hereafter MM96. These authors describe a model where the absorp- 
tion is associated with gaseous galactic halos in a two-phase structure: a shock heated hot 
phase and a photoionized phase of clouds moving through the halo. The authors associate 
this phenomenon with Lyman limit systems, though note that at high redshift [z > 2) the 
model predicts too few Lyman limits. An important aspect of this scenario is that strong 
Mgll absorption is primarily restricted to the halos of relatively massive galaxies. This is 
because the shock heated gas can cool only where the gas is dense enough to radiate effi- 
ciently, limiting the possible extent of the cold gas. Presumably the authors were motivated 
by observations showing that strong Mg II systems at z < 1 are frequently identified at im- 
pact parameters less than 30 hT 1 kpc from L pa L* galaxies (Lanzetta 1993; Steidel 1993). 
In any case, the MM96 scenario ties Mgll absorption to massive, passively evolving galaxies 
rather than merging or starburst events. 

This paper presents two lines of investigation into the nature of strong (equivalent width 1 
W r > 1 A) Mgll systems: (1) an automated search of the quasar spectra database of the 
Sloan Digital Sky Survey Data Release 3 (SDSS-DR3) for Mgll systems with z = 0.35 — 2.3; 
and (2) the analysis of a sample of Mgll systems (z = 0.8 — 2.7) identified in high resolution 
spectra of 92 quasars. The paper is organized as follows. § 2 explains the data used in this 
study and the analysis undertook on it. § 3 presents our results, including a plot of dN/dz 
versus redshift, a tabulation of each of the discovered systems, and a plot of interesting metal 
absorption lines. This section also includes a discussion of the relevance of our results to 
models of metal pollution of the IGM. The paper concludes in § 4 with a few remarks about 
future avenues of research. 



2. Data and Analysis 

This section describes two surveys for strong Mgll lines. In the first, the moderate 
resolution spectra of Data Release 3 from the Sloan Digital Sky Survey are searched for 
absorption. This survey provides the largest sample of Mgll systems compiled to date. The 
second data set is a much smaller sample of high resolution spectroscopy originally obtained 
to study high z damped Lya systems. This data extends the search to higher redshift and, 
more importantly, allows a quantitative assessment of the kinematics of Mg II absorbers. 



1 Unlcss otherwise stated, W r refers to the equivalent width of the Mgll A2796 transition. 
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2.1. SDSS-DR3 

The first set of data comes from Data Release 3 of the Sloan Digital Sky Survey (SDSS- 
DR3). The sample has been restricted to SDSS-DR3 quasars with redshift z qso > 0.35, 
for a total of 45023 quasars. SDSS spectroscopy has a resolution R m 2000 and covers 
3800 - 9200 A. This coverage implies a redshift search window of 0.35 < z < 2.3 for the 
Mgll AA2796,2803 doublet. 

A search and analysis of strong Mg II absorption was carried out in several steps. First 
each spectrum was continuum fit red-ward of Lyct emission using a Principal Component 
Analysis to fit the quasar emission lines and a b-spline algorithm to fit the underlying 
spectrum, which results roughly in a power-law with broad emission features. Second, all 
3.5a features were identified using a Gaussian filter method matched to the spectral resolution 
of the SDSS data. This method works by convolving a Gaussian with o = 1 pixels to both 
the flux and error arrays. The wavelengths where the convolved flux exceeds the convolved 
error array by 3.5 are considered significant absorption features and a line list is constructed 
for each quasar. Each resulting line list was searched for pairs of absorption features with 
separation matching the Mgll doublet. Finally, the equivalent width of the Mgll 2796 line 
was measured by directly summing the normalized flux. The equivalent width was measured 
by summing the difference between the continuum-fitted quasar flux and unity over a redshift 
centered box of width 3.589 A (half the spectral distance between the rest Mg 11(2796) and 
Mg 11(2803) lines). This final step is dominated by systematic uncertainty (in particular 
continuum fitting and contamination due to absorption by both systems at other redshifts 
and the sky), with an expected typical error of 0.2-0.3 A. To remove false-positives, a visual 
confirmation was then made for each system with measured W r > 0.8 A. 

To check the completeness of the survey a Monte-Carlo test was run on a representative 
sample of quasars. Artificial Mgll systems were added to quasar spectra, which were then 
searched for absorption as described above. The artificial systems were varied in strength 
and redshift and the quasars were selected at random from the complete SDSS-DR3 quasar 
sample. The results of this simulation are displayed in Figure 1. The solid line in the top 
panel reflect results for W r > 1.0 A systems versus redshift. The W r > 1.4 A results are 
displayed with a dashed line. The lower panel displays the W r of the systems as measured 
by the search algorithm versus the actual W r of the added system. The data have been 
binned and the values displayed are the average detected W r divided by the actual value, 
with the errorbars indicating the RMS of the systems in each bin. This figure indicates 
that this technique of absorption detection is > 95% complete for W r > 1.0 A systems, and 
provides fairly accurate measurements of W ri with typical errors of 0.1-0.2 A. A very small 
number of outliers in measured W r were detected in this analysis, primarily a result of sky- 
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line blending due to small errors in the detected system redshift. This effect was observed 
in very few systems, and is ignored in the following analysis. There is a slight systematic 
elevation of measured W r versus the added W r , which is an effect of contamination by real 
absorption systems in a given spectra and is ignored in the following analysis. To further 
assess the completeness and accuracy of the search algorithm, fifty random quasars were 
examined by eye. In this sample the automated search missed 2 out of 15 systems, each 
with W r < 0.7 A; the W r measurements of the automated search agreed well with this small 
sample of found systems. Because of Nyquist biasing, the number of discovered W r > 1.0 A 
systems is expected to be inflated hy W r < 1.0 A systems included due to error in the 
measurement of W r . As the phenomena considered in this work have no specific cut-off in 
W r , and as the effect is likely to be small, this error is neglected in the following analysis. 

In order to determine the number density of Mg II systems, one must calculate the total 
redshift path density g(z) of the survey. Large spectral regions, especially at higher redshifts, 
are highly contaminated by the sky. Simply removing these regions from the analysis would 
greatly reduce the redshift pathlength available to this study, but it is important to properly 
account for those regions which are too unreliable to accurately discover Mgll absorption. 
To account for this issue, the search was limited to regions where the Gaussian fitted (a 
= 1 pix) signal-to-noise (of the continuum) S/N > 7 per pixel, often resulting in non- 
contiguous search regions for a given quasar. The analysis of SS92 has been followed with 
a slight modification. For this work, instead of adding unity to a g(z) redshift bin for each 
spectrum which contains data for that bin, the decimal percentage of the given bin's redshift 
range actually covered in each spectra is added. Because strong Mg II systems dampen the 
noise of a spectra at their wavelength (i.e. Poisson noise is a significant source of error), and 
the S/N is calculated by dividing the continuum flux by the noise, the S/N at the specific 
wavelength of a given system will tend to increase over the local average. To protect against 
over-counting Mgll systems versus g{z), only systems discovered within regions in which 
the S/N smoothed over a 50 pixel box was greater than 7 are considered for the following 
analysis. A plot of g(z) for the SDSS search is presented in Figure 2. Limiting the search 
to smoothed spectral regions of S/N > 7 per pixel and a minimum W r threshold of 0.8 A, a 
total of 6511 Mgll systems were discovered, with 4835 of W r > 1.0 A, 2594 of W r > 1.4 A 
and 1329 of W r > 1.8 A. Another 3031 systems were discovered outside of S/N > 7 regions, 
2586 of which have W r > 1.0 A. 

A histogram of W r values for the SDSS-DR3 sample is plotted in Figure 3. The dis- 
tribution is steeply rising at W r = 1 A, inspiring confidence that the analysis is reasonably 
complete above this limit. A power- law fit to the W r > 1.0 A distribution, f{W r ) = CW' 5 , 
results in the parameters C = 490.4(1^) and 5 = 2.24(1^). This fit is plotted as the 
dashed line in Figure 3. The slope of this function is steeper than that measured in SS92 
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(5 = 1.65 ±0.03), but is more consistent with the measurements of others (Tytler et al. 1987, 
5 = 2.2 ±0.3). It is clear from the figure, however, that a power-law is not a good description 
of the observed distribution, under-predicting the number of systems at low W r and over- 
predicting the number of larger systems. The dotted line plotted in the figure represents a 
fit of a modified Schecter function to the W r distribution, of the form f(W r ) = BW^^e^ Wr , 
which is a much better fit to the data, particularly at high W r . The best fit values for 
this functional form are B = 2398(1^) and = 1.19 ± 0.07. It is likely that allowing the 
characteristic W r , W*, to be a fitted variable rather setting it to unity would improve the 
fit even further. Because this value has no bearing on the following discussion, this fitting 
and its interpretation are left to future work. All the fits were preformed using unbinned, 
Maximum Likelihood analysis. The errors are 95% confidence limits. 

Table 1 presents the SDSS quasar name, RA and DEC, emission redshift, the Mgll 
redshift, and W r of each W r > 1.0 A Mgll system discovered in the automated SDSS search, 
including those which passed the 'by-eye' confirmation but were eliminated from further 
analysis by the S/N > 7 per pixel criterion. 



2.2. High Resolution Observations 

The second set of quasar spectra used in this study includes 51 high S/N spectra 
(FWHM ~ 6.3-8.4 km s" 1 ) taken with the High Resolution Echelle Spectrometer (HIRES; 
Vogt et al. 1994) on the Keck I 10m telescope. These data were recorded and reduced as 
part of the UCSD/Keck I Damped Lya Abundance Database and a full discussion of the 
reduction process can be found in Prochaska et al. (2001). The high resolution sample also 
includes 41 R ~ 8000 spectra taken on the ESI spectrograph on Keck II. These data were 
collected as a part of the ESI/KECK II Damped Lja Abundance Database, details of which 
can be found in Prochaska et al. (2003). Table 2 summarizes the spectra searched in this 
study, presenting the name of each quasar, z em , data collection instrument (with 1 and 2 
referring to HIRES and ESI, respectively), the lower and upper limits of redshift available 
for Mg II detection, and information on discovered absorption systems. It should be noted 
that every spectrum in this data set contains a damped Lya system; these were ignored 
throughout the Mgll analysis. 

For each spectrum a 'by-eye' search was made for Mg II absorption systems. Given the 
strength of this doublet, these systems are conspicuous. After identifying a Mgll system, its 
redshift, equivalent width (Mgll 2796) and velocity width were measured, and the spectrum 
was searched for absorption due to other transitions associated with the system. The velocity 
width (Av) calculation followed the analysis of Prochaska and Wolfe (1997), i.e., the region 
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bounding 90% of the total optical depth in an unsaturated transition (typically Fe II 2600). 
The kinematic complexity of each system was also assessed, generally using complementary, 
non-saturated metal absorption lines (e.g. Fe II 2600). This measure is a simple summation of 
the distinct velocity components visible and was made due to the importance of absorption 
complexity to the arguments of BCCV. Because the ESI spectra are of lower resolution 
than those taken with HIRES, the observed complexity is expected to be reduced. Velocity 
structure will be lost for components separated by less than 45 km s -1 , especially in saturated 
absorption features. On the other hand, the spectra are sensitive to structure in the velocity 
profiles for components separated by greater than 50 km s _1 . 

The redshift path density g(z) has been calculated in a similar manner as the SDSS 
sightlines. The only significant difference is that spectral gaps inherent to the HIRES instru- 
ment add some uncertainty to this calculation. A possible prescription for accounting for 
the HIRES spectral gaps involves disqualifying those spectral areas immediately preceding 
and following the gaps, as only one half of the doublet is not enough to determine its nature. 
Because the systems of interest in this work are quite strong and because the wavelength 
coverage allows for investigation and confirmation of systems based on other metal absorp- 
tion lines (e.g. Fe II), this prescription was not followed and only the HIRES spectral gaps 
large enough to hide both components of a Mg II absorption doublet were disqualified in 
the calculation of g(z). Figure 4 presents g(z) for the high resolution sample. Clearly this 
high-resolution data extends to a higher redshift than available in the SDSS sample, though 
the small number of sightlines limits the significance of statistical results gleaned from this 
data. 

The high resolution survey is complete to below the W r = 1.0 A level. Both the ESI 
and HIRES datasets provide more than enough resolution and S/N for confidence in this 
statement. More than 50 Mgll systems with W r < 1 A were discovered in the search, though 
disregarded for the work presented here. At very long wavelengths (~ 10000 A in the ESI 
spectra), sky lines can confuse the search for Mgll systems. In a few particularly bad cases 
data are disqualified for this analysis, though the majority of spectra remained clear enough 
to easily identify absorption systems of the strength considered in this work. An example 
is presented in Figure 5: a plot of the Mgll system with W r = 0.56 A located at z — 2.51 
(~ 9800 A) in the spectrum of Q1337+11. Even crowded by the sky, this system is easily 
identified though it is only half as strong as the systems of interest to this study. 

The 22 W r > 1.0 A Mgll absorption systems discovered in this survey are depicted in 
Figure 6. For each system the Mgll 2796 absorption line is displayed in relative velocity 
(km s -1 ) space with km s" 1 centered on the redshift reported in Table 2. When more than 
one system was discovered in a given quasar spectrum, the lower and higher redshift systems 
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are labeled with a (1) and (2), respectively (no more than two systems were detected in any 
one spectrum). An unsaturated or mildly saturated transition (usually an Fe II line) is also 
displayed for each system. 

Six (~ 50%) of the systems for which spectra at the relevant wavelength were available 
(both observed in the quasar spectrum and not lost to the Lja forest) contained evidence of 
absorption due to Si IV, C IV, and Al III. These lines are displayed in Figure 7, where the 
systems are labeled as in Figure 6. These lines are displayed in relative velocity space with 
km s _1 centered at the system redshift found in Table 2. 

3. The Incidence of Strong Mg II Absorption 

3.1. Observations of and Fits to £ Mg (X) 

One means of investigating the physical nature of strong Mgll systems is to examine 
their incidence as a function of cosmological time. The most convenient and least ambiguous 
measure of this incidence is £ Mg (X), the line-density of absorbers per comoving Mpc. The 
incidence of a given class of absorbers at any redshift is proportional to their covering fraction 
on the sky. If the Mg II systems are a non-evolving population (i.e. constant number per 
comoving Mpc -3 n and constant cross-section a ), then one would predict ^Mg(A) = no<Jo. 
Should the number density or cross-section evolve with redshift, then no and a are replaced 
with n(z) and u{z). If, for example, strong Mgll systems are related to merger events, then 
one would expect n(z) to track the merger rate, predicted to be proportional to (1 + z) 3 (e.g., 
Patton et al. 1997), which implies a steep redshift dependence for £M g (X). Furthermore, 
because the number density evolution of halos in hierarchical cosmology is a sensitive function 
of mass, making comparisons with the evolution of £yt g (X) may indicate the mass of the 
systems giving rise to Mg II absorption. 

Historically, analysis has focused on the incidence of Mg II systems as a function of 
redshift, £Mg(z) (frequently written as dN/dz). Interpretation of this measure is somewhat 
ambiguous, because the evolution is influenced by the expansion of the universe as well as 
any evolution of the absorbers in question 

£ Mg ( z ) = £ Mg (X) * (c/H )(l + z) 2 [n M (l + zf + Ov]- 1 / 2 . (1) 

For this reason the majority of the analysis below is presented in terms of £ Mg (X) : which 
renders more obvious the trends discussed. For comparison to previous work, statistics of 
£Mg(z) are also presented. 
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Figure 8 presents ^Mg(A) for the Mgll systems identified in the SDSS-DR3 sample over 
the redshift range z = 0.35 — 2.3. The vertical errorbars assume Gaussian statistics while the 
horizontal errorbars indicate the size of each redshift bin. The SDSS results are displayed 
in three sets: the closed triangles represent the entire W r > 1.0 A dataset, the open squares 
indicate systems of 1.0 A < W r < 1.4 A, and the cross-hatches represent W r > 1.8 A. The 
data plotted in this figure are also listed in Table 3, as are the results for several other samples 
in W r . Qualitatively the samples shown in Figure 8 exhibit similar £ Mg (X) evolution. At 
redshifts z > 0.8, lyig(X) is nearly constant and linear fits yield slopes consistent with zero 
at the 95% c.l. At z < 0.8, in contrast, each sample exhibits a decrease in £ Mg (X). For the 
1.0 A < W r < 1.4 A sample the decline is marginally significant (~ la), while the effect is 
highly significant for any sample with W r > 1.4 A. 

We have considered two functional forms to describe the evolution of l Mg (X) with 
redshift: (1) a linear fit ^Mg(^) = b + mz and (2) an exponential fit ^M g (A) = iVexp(— z /z). 
Fitting a line to the l.oA < W r < 1.4A sample yields a slope consistent with zero: m = 
0.002±0.003. In contrast, the strongest Mg II absorbers (W r > 1.8A) have m = 0.011±0.002. 
Therefore, there is steeper increase in £yi g {X) for the larger W r absorbers, with the majority 
of the difference occurring at z < 1. These functional forms have very large reduced \ 2 
values, however, and are not good descriptions of the data. The exponential fits, in contrast, 
are a reasonable model of the data as these allow for a sharp dropoff of incidence at low 
redshift. Using standard maximum likelihood methods, we find N = 0.044 ± 0.003 and 
z = 0.12±0.07for the 1.0 A < W r < 1.4 A systems, N = 0.037±0.003 and z = 0.47±0.10 
for the W r > 1.8 A systems, and N = 0.11 ±0.006 and z = 0.28±0.05 for the full W r > 1.0 A 
dataset. These results are summarized in Table 4. 

For comparison to previous work, power-law fits to the ^Mg(^) values of the form N (l + 
z) 10 have also been performed. We stress, however, that such fits are no longer a good 
description of the data. In particular, this functional form cannot describe both the increase 
in ^MgPO at low redshift and the roughly constant value at z — 1 to 2. Nevertheless, this 
analysis was done for several cuts of the data. The maximum likelihood values for Nq and 70 
are N = 0.084, 7o = 1.40 for W r > 1.0 A, N = 0.051, 70 = 0.99 for l.oA W r > 1.4 A, and 
N = 0.016, 70 = 1.92 for W r > 1.8 A. Table 4 reports central values for N and 70, as well as 
95% confidence limits for each of the three fits. It should be noted that the values of 70 and N 
are highly correlated. Previous work on Mgll absorption has yielded similar, if somewhat 
steeper, results. SS92 find that for W r > 1.0 A systems, dN/dz oc (1 + 2 ) 2 - 24±0 - 76 . Also 
consistent with the results of SS92 is the increase in 7 with minimum equivalent width. We 
emphasize, however, that the differences in 70 are dominated by the evolution in £ Mg (X) at 
low redshift because all of the sub-samples have nearly constant ^Mg(A) at z > 1. Altogether, 
the W r > 1.0 A sample has an average incidence per unit redshift ^Mg(z) = 0.174 ± 0.002 



- 11 - 



for the redshift region 0.35 < z < 2.3. This value is consistent with but 25% lower than the 
value reported by SS92 for the redshift range 0.2 < z < 2.2 and are well matched to those 
of Nestor et al. (2004), who report results from the SDSS Early Data Release. 

3.2. Selection Bias 

It is important to understand and constrain possible biases in the SDSS-DR3 quasar 
sample. The value of £M g (X) is sensitive to observational bias in several forms. For example, 
if the Mgll systems are themselves massive enough to gravitationally lens, and thus brighten, 
their background quasars then the SDSS-DR3 would be statistically biased toward observing 
quasars which exhibit strong Mgll absorption (Menard 2005). This would lead to an 
underestimation of g(z) and an overestimation of £ug{X). The opposite effect may result 
if quasar spectra are reddened by dust in the Mg II host systems. Indeed, this bias is of 
particular concern, as many of the systems in this sample will be associated with damped 
Lya systems (Nestor et al. 2003) and such systems may cause reddening in the spectra of 
background quasars (Pei et al. 1991; Murphy & Liske 2004). 

Following the analysis of Murphy & Liske (2004) to address the matter of gravitational 
lensing, two averaged bootstrap samples of quasars free of strong Mgll absorption were 
created. The first of these samples was free of Mgll absorption in the range of 1.0 A 
< W r < 1.4 A, the second of W r > 1.8 A systems. These quasars where chosen to reproduce 
the redshift distribution of those sets of Mg II absorption-hosting quasars, but were otherwise 
chosen at random. Repeating this process many times uncovered the average non-absorption 
quasar magnitude distribution, which was then compared to the magnitude distribution of 
those quasars with Mg II absorption. This analysis revealed no signal of gravitational lensing 
in either the 1.0 A < W r < 1.4 A or the W r > 1.8 A datasets. Indeed, two-sided Kolmogorov- 
Smirnov tests reveal no evidence of differing magnitude distributions between the absorption 
and non-absorption datasets in either W r bin, with KS probabilities ranging from 0.2 — 0.6. 

Dust however, does seem to affect the spectra on average of those quasars hosting 
the strongest of the Mgll absorption systems. By building samples of quasars similar to 
the controls used in the lensing test (but now also constraining magnitude distributions of 
the control samples), composite 'non-absorption' quasar spectra were compiled, red- ward of 
Lya emission and in the rest wavelength of the quasar emission, for the two sets of W r Mg II 
systems described above. The second and fourth panels of Figure 9 display these spectra, as 
well as composite 'absorption' spectra, for the two Mgll absorption W r bins. Also plotted 
are power-law fits to these spectra, /(A) oc X a . Dust will preferentially absorb bluer light, 
tilting the spectrum and lowering the absolute value of a. 
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The first and third panels represent the result of dividing the composite Mgll spectra 
with the control spectra. The power-law fits to the 1.0 A < W r < 1.4 A composite spectrum 
and its companion control spectrum are nearly identical, |A(a)| ~ 0.02. These systems, 
therefore, are unlikely to impose an important dust bias. The larger W r sample clearly 
displays a dust signature, as seen in the divided spectrum in the upmost panel of Figure 9. For 
this sample |A(a)| ~ 0.11. Assuming Small Magellanic Cloud type extinction and employing 
the fitting formula of Pei (1992), this |A(a)| corresponds to a color excess E(B — V) ~ 0.01. 
This is an important result, especially as there are known quasar selection effects within 
SDSS, which uses colors in addition to apparent magnitude to select quasars for spectroscopic 
follow-up (Richards et al. 2002). Given the small color excess implied by the measured 
|A(a)|, however, these effects are likely to small. A full consideration of these effects will 
be the subject of a future paper. For now, we note that the effect of reddening would 
generally be to under-measure ^Mg(A). The net effect could be a systematic shallowing of 
^Mg(A) especially at z > 1, although the reddening is sufficiently small that the corrections 
to ^Mg(A) should be minor. 



3.3. Implications and Interpretations 

The evolution of ^Mg(A) described by Figure 8 has several implications for the physical 
processes relevant to Mg II systems. First, the near constancy of ^Mg(A) over the ~ 3 Gyr 
spanning z = 2 to 0.8 indicates the cross-section to strong Mg II systems is very likely to be 
dominated by dark matter halos with masses M <C 1O 12 M . Although structure forms much 
earlier in the ACDM cosmology than an Einstein de Sitter model, the number density of 
galaxies with M > M*(z = 2) « 1O 11 M increases exponentially over this epoch (e.g. Jenkins 
et al. 2001). If these massive galaxies dominated the cross-section, then a would have to 
decrease exponentially to yield constant £ug(X). We consider such a coincidence to be highly 
unlikely. It is far more reasonable to physically associate the phenomenon with galaxies 
having mass M < M*(z = 2) and assume a nearly constant cross-section during this epoch. 
Second, the decline in ^Mg(A) at z < 0.8 indicates the cross-section to strong Mg II systems 
is decreasing in time. This conclusion applies particularly to the W r > 1.8 A absorbers where 
the decrease in ^Mg(^) is secure. For the same argument as above, the number density of 
dark matter halos is either constant (for M <C M*) or increasing substantially (for M 3> M*). 
A decrease in £ Mg (X) therefore, can only be explained by a decrease in a. The observations 
imply the processes responsible for strong Mg II absorbers are 'turning off' at z ~ 1. Third, 
one identifies a rough correspondence between the star formation history of the universe 
(Madau et al. 1996) and ^Mg(A). Granted the two quantities track one another over the 
redshift range z = 0.5 to 2, it is tempting to connect the strong Mg II systems with processes 



-13- 



related to active star formation (e.g. galactic winds, starbursts). 

We now examine the implications of the evolution of ^Mg(A) hi terms of the two models 
described in the introduction: (1) a starburst scenario where the incidence and kinematics 
of strong Mg II systems are determined by feedback processes (BCVS); and (2) infalling 
gas in the outer halos of massive galaxies (MM96). An important prediction of the MM96 
scenario is that £Mg(X) should steeply decline at high redshift once the characteristic dark 
matter halo mass M* becomes smaller than the typical halo mass giving rise to strong Mg II 
absorption. In the MM96 scenario, the halo mass is ~ 1O 12 M , which is necessary to insure 
the halo has a significant reservoir of cool, photoionized gas (MM96). Within the ACDM 
cosmology with <t 8 = 0.9, M* = 10 12 M Q at z = 0.75. Therefore, this model predicts a steep 
decline in £ Mg (X) at z > 1 which is not observed. As noted above, to adopt a model where 
strong Mg II absorption arises in an unevolving population of dark matter halos, the typical 
mass would need to be less than M*(z — 2) « lO ia5 M . Therefore, we conclude that the 
majority of strong Mg II absorbers do not arise in massive halos typical of the W r ~ 0.3A 
population at z ~ 1. In fairness, we note that the MM96 model was primarily introduced 
to explain the W r ~ 0.3A Mg II systems which may have a different origin than the strong 
Mg II absorbers considered in this paper. 

It is more difficult to compare the observations against the starburst scenario because 
there does not exist a well-developed model. Nevertheless, consider the following arguments. 
At z ~ 1, the dark matter halo merger rate is expected to scale as (1 + z) 3 (Patton et al. 
1997). If one assumes the starburst phenomenon is merger- driven, then the simple prediction 
is an increase in £ivig(A) with increasing redshift. Indeed, this is observed from z = 0.3 to 0.8, 
yet all of the sub-samples are inconsistent with even a (1 + z) increase in ^Mg(A) at z > 1. 
Because the incidence is the product of the number density and cross-section of the halos, 
one could envision a scenario where a is decreasing with redshift as n increases such that 
£ Mg (X) remains roughly constant. For example, the Mg abundance could be much lower 
at higher redshift or the velocity fields more quiescent, but these explanations are rather 
unlikely for a starburst model. Recent observational results of Len et al. (2004) suggest 
another possible solution. These authors find that the galaxy merger rate evolves much 
more shallowly with redshift than dark matter haloes, (1 + /2 )°- 51±0 - 28 ; a much better match 
to the Mg II observations. 

As noted above, the general trend of ^Mg(A) roughly tracks the global density of star 
formation. That is, the star formation rate (SFR) density is described as roughly constant 
from z — 1 to 2 (where it is not well constrained) and decreases by a factor of several from 
z — 1 to today. Although the correspondence between the incidence of strong Mg II absorbers 
and the SFR density could be simple coincidence, it may be expected if the absorption 
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systems arise from feedback processes related to star formation. The correspondence may 
suggest that significant star formation is required to maintain a large cross-section to strong 
Mg II absorption, i.e. this phase is short-lived in galactic halos and requires continued input 
to maintain the observed incidence at z — 1 to 2. At the very least, we consider the similarity 
to be at least anecdotal evidence in favor of a starburst scenario. 

The actual value of the incidence of strong Mg II absorption - as opposed to its evolution 
- provides an important constraint on the hypotheses of massive galactic halos and starbursts. 
The phenomenology of MM96 implies a typical impact parameter for Mg II absorption of 
p ~ 30 h -1 kpc. Using ^Mg(AT) ~ 3 x 10~ 5 h Mpc" 1 at z ~ 1, the value measured in 
this work for W r > 1.0 A systems, a number density n ~ 0.01 h 3 Mpc 3 is measured for 
large, quiescent galaxies at this redshift. It should be recalled that the measured £ Mg (X) for 
Mg II at this redshift may be somewhat influenced by merging systems, so the value used for 
£Mg(X) in this calculation should be considered an upper-limit. Again assuming standard 
concordance cosmology and using the mass function introduced by Jenkins et al. (2001), 
the mass of halos with n ~ 0.01 — 0.03 h 3 Mpc -3 is M ~ 10 11 M . Similar to our previous 
conclusions, we find that the MM96 scenario is nonviable at z > 1 if one adopts their typical 
mass halo of 1O 12 M unless one extends the impact parameter to an unlikely p ~ 100 kpc. 

Turning to the starburst hypothesis of BCCV, the observations can be used with pre- 
vious work on superwinds and LBGs to constrain the required impact parameter of the 
outflows. If high redshift Mgll systems do correspond to SB galaxies, one expects them 
to trace LBGs, also thought to harbor starbursts. Using the measured number density 
n ~ 0.016 h 3 Mpc -3 for LBGs at z ~ 3 (Steidel et al. 1999), and assuming the value of 
^M g (^0 observed at z = 2 for the W r > 1 A sample (2 x 10~ 5 h Mpc -1 ), the predicted o for 
merging Mgll systems at z = 3 is a ~ 10~ 3 hr 2 Mpc 2 . The value of £ Mg (X) used in this 
calculation may be considered a maximum, as the overall sample of Mg II systems discovered 
in this survey is likely to be polluted with large-scale processes, Damped Lja absorption 
systems, or other systems, and thus the resultant value of a is also a maximum. To determine 
whether this value is consistent with the wind hypothesis, theoretical upper and lower bounds 
on a may be used. Assuming a tw^e ~ 100 km s" 1 and a bubble age t r ~ 60 Myr (BCCV, 
a lower limit equal to the expected time-scale of the SB), the approximate impact parameter 
will be a ~ 1 x 10~ 4 Mpc 2 . The expansion time and velocity used for this calculation are 
both lower limits, so this is a lower limit. Computer simulations presented in Furlanetto & 
Loeb (2003) find the maximum size of these bubbles (after very long evolutionary periods) 
to correspond to a ~ 3 x 10~ 2 Mpc 2 . The gas in such large bubbles, however, may be of too 
low density to give rise to strong Mg II absorption, so this should be considered an extreme 
upper limit. In any case, the cross-section determined from Mgll absorption and LBG data 
is consistent with the theory that superwinds are the primary hosts of high redshift Mgll 
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absorption. 

4. Kinematics 

Examining the profiles in Figure 6, it is evident that the majority demonstrate at least 
one of the three criteria discussed in BCCV for the identification of superwind driven galactic 
bubbles. Because these features were clearly present in even the weaker W r systems, the 
BCCV criterion of W r > 1.8 A has been relaxed to W r > 1.0 A. Based primarily on 
unsaturated transitions (usually a Fell line), the majority (~ 80%) of the discovered systems 
exhibit some degree of complexity, though as discussed above, kinematic features of width 
Av < 45 km s _1 are unobservable in the ESI spectra. Because every system discovered in a 
HIRES spectrum displays a great deal of spectral complexity, and as many of these features 
are kinematically closer than the above specified limit, it is probably a safe assumption that 
most or all of the systems studied in this work conform to the complexity criterion introduced 
out by BCCV. 

Figure 10 presents a histogram of the velocity widths of the high resolution sample. 
Clear is the cutoff at ~ 100 km s _1 which corresponds to the W r > 1.0 A imposed on this 
study (although recall that the At> measurements are based on unsaturated profiles). The 
velocities range from Av = 66 km s -1 for the z = 1.328 system toward Q1331+17 to Av = 
470 km s _1 for the z = 1.339 system toward PSS2241+1352, with the majority of systems 
having Av ~ 200 km s -1 . In the idealized case of a spherical shell expanding with constant 
velocity v s h e ii, the velocity width for an impact parameter of R/2 is \^3v s h e u- Therefore, the 
median velocity width for this sample corresponds to a shell speed of v s h e ii ~ 150 km s _1 
. This value is consistent, though on the low end, with numbers observed for local galactic 
superwinds (Heckman et al. 2000) and those predicted from numerical simulations (e.g., 
Strickland & Stevens 2000). An important and somewhat unexpected result of this analysis 
is that very few Mgll systems show Av > 250 km s" 1 . If this gas is tracing winds arising 
from starburst systems, it is evident that the optical depth of rapidly expanding 'bubbles' 
at z ~ 1.5 is small: r < 1.0 x 10~ 5 h Mpc" 1 (95% c.L). 

In one respect, the data displayed in Figure 7 are generally unsupportive of the Mg 
II absorber-superwind hypothesis. In a wind model, one envisions a cold, expanding shell 
of gas which envelopes a rarefied, hot region which in turn surrounds the post-SB galaxy. 
In this scenario, one would not expect to observe high-ionization state particles (e.g. Si 
IV) in the same distribution as the colder gas (Mg II; Strickland & Stevens (2000)). One 
might expect the higher-ionization state gas to be distributed in a kinematically narrower 
region (it is no longer expanding). Alternatively, the gas may exhibit turbulent velocities 
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at the virial velocity dispersion of the galaxy and not conform to the clump patterns of 
the gas in the shell. In either case, these trends are not observed in those systems with 
high-ionization lines. Figure 7 shows a clear trend for high- and low-ionization state gas to 
share the same regions of km s _1 space and to often display similar clumping patterns. This 
result is consistent with the work of Churchill et al. (1999), who reported strong kinematic 
links between Mgll and CIV absorption in the redshift range 0.4 < z < 1.4, and discussed 
consequences for both the above scenario and that proposed by MM96, which involves cold 
bubbles of gas confined by pressure in a hot gas medium. 

The construction of MM96 also implies different kinematic situations for the hot and 
cold gas, a prediction not supported by the data. Indeed, the models presented in MM96 
predict that systems showing strong Mg II absorption (large galaxies) should contain little if 
any C IV absorption, which is clearly ruled out by the data. A suggestion put forward to fix 
this discrepancy is a "warm" phase of collisionally ionized gas in an interface between the 
cold bubbles and the hot medium. This would explain the observed kinematic link between 
Mg II and more highly ionized species (in either the expanding bubble or two-phase gaseous 
halo scenarios), though the authors point out that the amount of collisionally ionized CIV 
found in such a phase would likely be small. There is also the possibility that CIV arises 
at higher radii than the simple models of MM96 consider. This solution, while providing for 
more C IV absorption, does not provide a simple explanation for the observed kinematic link 
between species. More sophisticated models of the two-phase halo structure are needed to 
clarify this issue. 



5. Concluding Remarks 

There are many open questions left in the task of assigning a particular phenomenology 
to strong Mg II absorption. Detailed simulations of superwinds would be useful in pinning 
down the expected absorption profile of such a system, as well as providing information 
on expected metallicities and temperatures. A handful of strong Mg II absorption-selected 
galaxies have been studied kinematically (see Steidel et al. 2002; Churchill et al. 2000), 
with varying results, though the number remains small and little consideration has gone into 
the link stipulated in this work. A more comprehensive survey over a larger redshift range 
would allow for stronger statements regarding the correlation between these absorbers and 
galaxy populations, as well as a refined understanding of their kinematic nature. Optical 
studies could only provide results at low redshift, but confirmation of the SB and post- 
SB nature expected of the hosting galaxies, even for low-z examples, would strengthen the 
argument considerably. Adaptive optics might be useful here, and a program has been begun 
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at Lick Observatory to obtain AO images of a few quasars known to harbor low redshift, 
strong Mg II systems in an attempt to morphologically identify these systems. We are also 
performing a separate imaging program with ACS on the Hubble Space Telescope (Cycle 14) 
to cover a wider range of redshift. Finally, with the advent of infrared spectrographs (e.g. 
the GNIRS instrument on Gemini South) and the availability of large numbers of bright, 
high redshift quasars (as in the SDSS DR3), the phenomenon can be studied to very high 
redshift (z ~ 4). In particular, it is critical to further test the assertion that strong Mg II 
systems trace processes related to star formation. 
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Fig. 1. — Completeness measurements for the automated SDSS-DR3 absorption search, 
based on Monte-Carlo simulations (S/N > 7). The top panel indicates percent completeness 
versus absorption redshift, with the solid line indicating systems of W r > 1.0 A and the 
dashed line for systems of W r > 1.4 A. The bottom panel depicts the W r measured by the 
automated search algorithm versus the actual W r of the system added, with the vertical axis 
indicating the average ratio plus RMS error of systems for each 0.1 A W r bin (the horizontal 
axis) . 
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Fig. 2. — The redshift path density, g{z), probed by SDSS-DR3. This represents the total 
number of lines of sight at a given redshift in which it is possible to detect a strong Mg II 
system in the SDSS-DR3 spectra. Clear in this plot are sky-absorption features at high 
redshift, as well as very strong lines at z ~ 1 at z ~ 1.25. 
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Fig. 3. — A histogram of the W r values for the Mgll systems discovered in SDSS-DR3 
with our automated search algorithm. The vertical dashed line is located at W r = 1.0 A. 
Given that the incidence of W r continues to rise well below lA, the survey is expected to 
be > 95% complete. A power law fit for systems of W r > l.oA is plotted as a dashed line; 
f(W r ) = 490. 37W 7 "" 2 - 245 . The dotted line represents a modified Schecter functional fit the 
the W r > 1.0A data, f(W r ) = 2398.4W- L19 e -^. 
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Fig. 4. — The redshift path density, g(z), probed in the high resolution survey. The dashed 
vertical lines indicate the redshift region considered in the work presented in this paper. 
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Fig. 5. — An example Mgll system discovered at z — 2.51 in the spectrum of Q1337+11. 
This system is located at ~ 9800A, where the sky begins to crowd the spectrum yet is easily 
identified at this spectral resolution and signal-to-noise. 
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Fig. 6. — The 22 detected W r > 1.0 A, Mgll absorption systems detected in the high 
resolution survey. For each system the Mgll 2796 absorption line is displayed in relative 
velocity space with km s _1 centered at the redshift reported in Table 2, with (1) referring 
to the low-z and (2) to the high-z system in cases of multiple detections per quasar. For 
each system, an associated, less saturated metal transition is displayed below the Mgll 2796 
feature. The systems are ordered by increasing redshift from top to bottom then left to 
right. 
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Fig. 7. — Six systems for which transitions of Al III, C IV, and Si IV are detected. This figure 
adopts the same labeling as Figure 6. 
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Fig. 8. — The incidence of strong Mg II absorbers, £ug(X), from the SDSS-DR3 over the 
range 0.35 < z < 2. The triangles show the full sample (W r > lA), the open squares are all 
systems of 1.0 A < W r < 1.4 A and systems of W r > 1.8 A are indicated with cross-hatches. 
The differences between the datasets in this plot have important implications for the nature 
of Mg II absorption hosts. 
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Fig. 9. — Results from an analysis of the dust reddening of quasar spectra hosting strong 
Mgll absorption. The top two panels are results from the W r > 1.8 A systems with the 
lower panel displaying composite spectra of quasars hosting strong Mg II absorption (black) 
and quasars without absorption (red). The top panel is the absorption spectra divided by 
the non-absorption spectra. The bottom two panels are the equivalent results for systems of 
1.0 A < W r < 1.4 A. These results are indicative of the spectra of those quasars hosting the 
largest of Mg II systems having been reddened. 
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Fig. 10. — A histogram of Av for the 22 systems discussed in this paper. Note the sharp 
cutoff at ~ 100 km s _1 corresponds to the W r > 1.0 A limit imposed for consideration of a 
Mgll system. More importantly note the near absence of systems with very large Av value. 



-33- 



Table 1. SDSS-DR3 W r (279Q) > 1.0 A Mgll Systems 



Quasar 


RA 


Dec 


Zqso 


ZMg II 


W r (2796) 


J000009.42-102751.87 


00:00:09.42 


-10:27:51.87 


1.844 


1.313 


1.69 


J000050.60-102155.87 


00:00:50.60 


-10:21:55.87 


2.640 


1.056 


1.17 


J000053. 72+150637.20 


00:00:53.72 


+ 15:06:37.20 


1.411 


0.897 


1.63 


J000143.41+152021.49 


00:01:43.41 


+ 15:20:21.49 


2.638 


1.598 


1.03 


J000221. 11+002149.46 


00:02:21.11 


+00:21:49.46 


3.057 


1.958 


2.33 


J000221.80+151454.58 


00:02:21.80 


+15:14:54.58 


1.823 


1.172 


1.05 


J000221.80+151454.58 


00:02:21.80 


+15:14:54.58 


1.823 


1.435 


1.29 


J000221.80+151454.58 


00:02:21.80 


+15:14:54.58 


1.823 


0.716 


1.22 


J000221.80+151454.58 


00:02:21.80 


+15:14:54.58 


1.823 


0.846 


1.68 


J000341. 22-094101. 03 


00:03:41.22 


-09:41:01.03 


1.739 


1.584 


2.39 



Note. - The complete version of this table is in the electronic edition of the 
Journal. The printed edition contains only a sample. 
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Table 2. High Resolution Survey Quasars and Found Mg II Systems 



QSO 


Inst 








nin 




nax 


ZMgll 


Av 


W r Metal; 


s components 


RDflfll Q 1 K 

.orvuuiy-io 


i 
i 


A 
4. 


^3 
JO 


1 

1 . 


4D 

4U 


1 

1 


.yo 










RRDQ^I D4 

UIXUc/O -L U4 


i 


A 

4 


37 

o t 


1 

1 . 


3fi 
. OU 



z 


nn 

.uu 










uqn7A1 _i_47 


i 
i 


Q 

o. 


99 
zz 


u. 


81 


1 
1 


fi7 
.U ( 










109^4-00 


i 
i 


Q 
O. 


Q7 
. y i 


i 

i . 


1 fi 

. 1U 


1 

1 


Q9 
.yz 










PNQ^7 
rnyo / 


1 9 


o 
z 


uy 


n 

u. 


47 


9 
z 


7fi 
. / u 










roouyo i too 


1 9 


A 
4. 


9^ 

zo 


1 

1 . 


9Q 

zy 


9 
Z 


.00 












1 
l 


4:. 


41 


1 

1 . 


38 
oo 


9 
z 


1 9 
. 1 z 










nnnnn 9fi 
v^>uuuu-zu 


1 
1 


/I 


1 1 
. 1 1 


1 

1 . 


9fi 
ZU 


i 
1 


Q1 
.y i 












1 
1 


9 
z. 


.40 


U. 


4fi 

4U 


i 
1 


33 
.oo 










nn9m 4-1 1 


1 
l 


Q 
O. 


fi1 


i 

i . 


Ul 


1 
i 


73 
. i o 












1 

_L 


9 
z 


4Q 


n 

u . 


U J 


1 


.o / 










00940 99 
v^uz4y-zz 


1 
1 


Q 
O 


9n 
zu 


n 
u 


3n 
ou 


n 
u 


.OU 










nn^99 ^9 

Vc^UOZZ-OZ 


1 
1 






n 

u. 


49 

. 4Z 


n 

u 


Q9 
.yz 












1 
1 


Q 
O 


on 
.zu 


n 

u. 


.ou 


1 

1 


98 

.Zo 










^UOI ( -OO 


1 
1 


Q 
O. 


93 
zo 


n 

U. 


84 


1 

1 


R9 
.uz 


1 4^7 
1 .40 ( 


91 (\ 
Z 1U 


9 937 
Z.ZO ( 


7 


nnA^n 1 3 


1 
1 


9 
z. 


3n 

ou 


u. 


34 

.04 


n 

u 


7fi 
. ( u 










nn/i n9 

^U40o-UZ 


1 
1 


9 
z. 


9Q 

zy 


n 
u. 


/I 1 
41 


1 
1 


9^ 
.zo 










^uoor-oo 


1 
1 


9 
z. 


39 

oz 


n 
u. 


49 
.4Z 


1 
1 


93 
.Zo 










nn7 c i7-i- c i9 


1 
1 






n 
u. 


39 
.oz 


n 

u 


83 

.oo 










ons3fi4-i 1 


1 
1 


9 
z. 


7n 


n 
u. 


^9 
oz 


1 

1 


^7 

.0 ( 










ons4i 4-1 9 


1 
1 


9 
z 


90 
zu 


n 

u. 


71 


1 

1 


^7 

.O i 


1 DQ8 


9DR 
zuo 


1 4^0 

1.40U 


4 


^Uo41-|-lZ 
















1 91 Q 

i.ziy 


97/1 
Z I 4 
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Q0930+28 
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3. 


.42 


0. 


.86 


1 


.56 










Q0952-01 


1 


4. 


.43 


1. 


36 


1 


.91 










Q1005+36 


1 






0. 


.80 


1 


.51 










Q1055+46 


1 


4. 


13 


1. 


.25 


1 


.86 


1.386 


210 


1.164 


4 


Q1104-18 


1 


2. 


31 


0. 


.44 


1 


.47 










Q1108-07 


1 


3. 


.92 


1. 


.12 


1 


.98 










Q1202-07 


1 


4, 


69 


1. 


.49 


2 


.00 










Q1210+17 


1 


2, 


.54 


0. 


.56 


1 


.21 










Q1215+33 


1 


2. 


.61 


0. 


39 


1 


.26 
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Table 2 — Continued 



QSO 


Inst 


3 


QSO 




TL1TI 


Zmax 


ZMg II 


Av 


W r 


Metals components 


Q1223+17 


1 


2. 


92 





.71 


1 


.57 














Q1331+17 


1 


2. 


.08 





.51 


1 


.37 


1 


.328 


66 


1.189 




2 


Q1346-03 


1 


3. 


99 


1 


.18 


1 


.98 














Q1425+60 


1 


3. 


.17 





.81 


1 


.21 














Q1759+75 


1 


3. 


.05 


1 


.14 


1 


.71 














Q1850+40 


1 


2. 


.12 


0. 


.82 


1 


.79 














Q2038-01 


1 


2. 


.78 





.65 


1 


.21 














Q2206-19 


1 


2. 


.56 





.52 


1 


.33 














Q2223+20 


1,2 


3. 


,56 





.99 


2 


.65 














Q2230+02 


1 


2. 


.51 





.38 


1 


.46 


1 


.059 


110 


1.234 




4 


Q2231-00 


1 


3. 


02 





.74 


1 


.29 














Q2233+13 


1 


3. 


.27 





.86 


1 


.61 














Q2343+12 


1 


2. 


.52 


0. 


.43 


1 


.28 














Q2344+12 


1 


2. 


.79 





.55 


1 


.17 














Q2348-01 


1 


3. 


.01 


0. 


.81 


1 


.68 














Q2358-02 


1 









,66 


1 


.33 














BQ1021+30 


2 


3. 


.12 


0. 


.73 


2 


.65 


2 


.329 


200 


1.023 


SiIV,CIV,AlIII 


3 


BRI1 144-07 


2 


4. 


.16 


1 


.26 


2 


.65 


1 


.908 


190 


1.833 




2 


BRJ0426-2202 


2 


4. 


30 


1 


.32 


2 


.65 














CTQ460 


2 


3. 


.13 





.68 


2 


.65 














FJ0747+2739 


2 


4. 


.11 


1 


.25 


2 


.65 














FJ0812+32 


2 


2. 


.70 





.63 


2 


.65 


1 


.223 


260 


1.200 




3 


FJ0812+32 
















2. 


.626 


280 


1.300 


SiIV,CIV,AlIII 


1 


FJ2129+00 


2 


2. 


.94 


0. 


.72 


2 


.65 














FJ2334-09 


2 


3. 


.33 





.90 


2 


.65 


2 


.152 


140 


1.363 


SiIV,CIV,AlIII 


1 


HS1132+22 


2 


2. 


.88 


0. 


.66 


2 


.65 














HS1437+30 


2 


2. 


99 





.70 


2 


.65 














PC0953+47 


2 


4. 


.46 


1 


.40 


2 


.50 














PKS1354-17 


2 


3. 


.15 


0. 


.79 


2 


.65 














PSS0106+2601 


2 


4. 


32 


1 


.32 


2 


.65 


2. 


.198 


130 


1.110 




1 


PSS0133+0400 


2 


4. 


.13 


1 


.25 


2 


.65 


1 


.666 


160 


1.632 




2 
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Table 2 — Continued 
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X . Zj(J I 
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Zj 
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9 
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1 ^9 
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9 fi5 
z.uo 


1 33Q 
i .ooy 
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^ 1 u 


4 R9fi 




o 
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PQQ9941 4-1 359 










9 549 

Z.O^iZ 


1 70 
1 1 u 


1 779 

1. MZ 




9 
z 


PSS2323+2758 


2 


4.18 


1.25 


2.65 












PSS2344+0342 


2 


4.30 


1.29 


2.65 












PX0034+16 


2 


4.29 


1.32 


2.65 


1.800 


160 


1.576 




2 


Q0112-30 


2 


2.96 


0.71 


2.65 












Q0821+31 


2 


2.61 


0.61 


2.65 


2.535 


190 


1.738 


SiIV,CIV,AlIII 


2 


Q0933+28 


2 




0.93 


2.65 












Q1209+09 


2 


3.30 


0.86 


2.65 


2.585 


270 


3.020 


SiIV,CIV,AlIII 


2 


Q1337+11 


2 


2.92 


0.68 


2.65 












Q1502+48 


2 


3.20 


0.82 


2.65 


2.272 


220 


1.807 


A1III 


2 


Q2342+34 


2 


3.01 


0.77 


2.65 












SDSS0127-00 


2 


4.06 


1.22 


2.65 












SDSS2350-00 


2 


3.01 


0.75 


2.65 


0.863 


200 


1.474 




3 


SDSS2350-00 










2.425 


280 


2.732 


SiIV,CIV,AlIII 


2 
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Table3. SDSS-DR3 Mgll dN/dX 



<z> dN/dX N dN/dX N dN/dX N 

w r > 1.0 A i.oA<^< 1.4 A w 7 ,, > 1.8 A 



0. 


.41 


0. 


.056 


± 


0, 


006 


96 





.030 


± 





.004 


52 





.009 ± 


0, 


.002 


16 


0. 


.50 


0, 


,059 


± 


0, 


004 


190 





.030 


± 





.003 


97 





,012 ± 





.002 


41 


0. 


,60 


0. 


066 


± 


0, 


004 


253 





.033 


± 





.003 


128 





.016 ± 





,002 


61 





,70 


0, 


,081 


± 


0. 


004 


351 


0, 


.036 


± 


0, 


.003 


157 





,025 ± 





,002 


109 





.79 


0. 


,088 


1 

± 


0, 


004 


403 





.045 


1 

± 





.003 


206 





.021 ± 





,002 


96 





,89 


0, 


,089 


± 


0, 


004 


417 


0, 


.045 


± 





.003 


211 





,023 ± 





,002 


106 





.98 


0. 


,080 


± 


0, 


004 


328 





.036 


± 





.003 


147 





.022 ± 





.002 


91 


1 


.08 


0. 


091 


± 


0. 


,005 


297 


0, 


.038 


± 


0, 


.003 


124 





,028 ± 





,003 


91 


1 


.18 


0, 


092 


± 


0, 


005 


396 





.043 


± 





.003 


187 





,026 ± 





,002 


112 


1 


.27 


0. 


,089 


± 


0, 


004 


388 





.043 


± 





.003 


185 


0, 


.023 ± 


0, 


.002 


101 


1 


,37 


0. 


,102 


± 


0. 


,005 


451 


0. 


.048 


± 


0, 


.003 


215 





,026 ± 





,002 


116 


1 


.46 


0. 


,088 


± 


0. 


005 


336 





.039 


± 





.003 


149 





.023 ± 


0, 


.002 


89 


1 


,56 


0. 


,097 


± 


0, 


005 


315 





.040 


± 





.003 


130 





.034 ± 





,003 


110 


1 


,66 


0. 


100 


± 


0. 


006 


249 





.044 


± 





.004 


108 





,030 ± 





,003 


74 


1 


,75 


0, 


,084 


± 


0. 


,007 


129 


0, 


.036 


± 


0, 


.005 


55 





,023 ± 





,004 


35 


1 


,85 


0. 


,089 


± 


0, 


009 


100 





.039 


± 





.006 


44 





,027 ± 





,005 


30 


1 


,94 


0. 


,094 


± 


0. 


010 


84 


0, 


.036 


± 





.006 


32 





.035 ± 





,006 


31 


2 


.04 


0. 


,076 


± 


0, 


013 


33 





.016 


± 





.006 


7 





.032 ± 





.009 


14 


2. 


.14 


0. 


,070 


± 


0, 


020 


12 





.012 


± 





.008 


2 


0, 


.029 ± 


0, 


.013 


5 


2. 


.23 


0, 


065 


± 


0, 


,024 


7 





.046 


± 





.021 


5 





,009 ± 





.009 


1 



Note. — All dN/dX values reported are multiplied by c/Hq for ease of 
readability, and are in units of h Mpc^ 1 . 
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Table 4. dN/dX and dN/dz Functional Fits to SDSS-DR3 Mg II Systems 



W r (k) 


Shape" 


Normalization 




dN/dX = N + mz 


> l.U 


0.022 ± 0.004 


0.036 ±0.004 


> 1.0 


0.011 ±0.002 


0.011 ±0.002 


l.U < W r < 1.4 


0.002 ±0.003 


0.036 ±0.002 




dN/dX = iVexp" 


-zo/z 


> 1.0 


0.275 ±0.05 


0.113 ±0.006 


> 1.8 


0.467 ±0.10 


0.037 ±0.003 


1 < W,. < 1 4 


0.123 ±0.07 


0.044 ± 0.003 




dN/dX = N(l + 


z y< 


> 1.0 


0.54 ±0.10 


0.06 ±0.004 


> 1.8 


0.98 ±0.154 


0.01 ±0.001 


1.0 < W r < 1.4 


0.13 ±0.15 


0.04 ±0.004 




dN/dz = N(l + 




> 1.0 


1.40(1.24,1.56) 


0.08(0.095,0.075) 


> 1.4 


1.74(1.52,1.96) 


0.036(0.042,0.030) 


> 1.8 


1.92(1.60,2.22) 


0.016(0.021,0.013) 


1.0 < W r < 1.4 


0.99(0.77,1.28) 


0.051(0.060,0.041) 


1.4 < W r < 1.8 


1.56(1.25,1.89) 


0.020(0.025,0.015) 


m, zq, or 7 



Note. - - These are the la errors for the minimized- 
X 2 power-law fits and the 95% Confidence Limits for the 
maximum likelihood power-law fits to the SDSS-DR3 
dN/dz data. For the dN/dz results, the upper and lower 
Confidence Limits for 7 and iV are presented in paren- 
theses. The values of iV for these fits are calculated by 
constraining the integral of g(z) * dN/dz to equal the 
observed number of systems for a given cut in W r . 



